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We have studied the effects of a disordered optical potential on the transport and phase coherence 
of a Bose-Einstein condensate (BEC) of ^Li atoms. At moderate disorder strengths (Vd), we ob- 
serve inhibited transport and damping of dipole excitations, while in time-of- flight images, random 
but reproducible interference patterns are observed. In- situ images reveal that the appearance of 
interference is correlated with density modulation, without complete fragmentation. At higher Vd, 
the interference contrast diminishes as the BEC fragments into multiple pieces with little phase 
coherence. 
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The behavior of a superfluid/superconductor in the 
presence of disorder is of fundamental interest. A su- 
perfluid can flow without friction around obstacles and a 
superconductor can have zero resistance despite material 
defects. On the other hand, disorder is able to local- 
ize particles, resulting in an insulating state pL^. Experi- 
mentally, disorder-induced superfluid/superconductor to 
insulator transitions (SIT) have been probed in many 
systems, including superfluid helium in porous media [2], 
thin film and granular superconductors [H and ran- 
dom Josephson junction arrays [5]. While many believe 
that such a SIT is a quantum phase transition driven 
by quantum fluctuations, it remains a central task to 
understand exactly how the superfluid/superconducting 
order parameter, which consists both an amplitude and 
a phase, may be destroyed with increasing disorder. Nu- 
merous fundamental questions remain, such as the nature 
of the insulator, the fate of phase coherence throughout 
the transition, and the possibility of intermediate metal- 
lic phases [3l El [7] . 

Cold atoms, with their intrinsic cleanliness coupled 
with remarkable controllability of physical parameters, 
have emerged as exceptional systems to study various 
condensed matter problems. Recently, several experi- 
ments [8] m [TOl [11] have studied ^^Rb condensates in 
random optical potentials, and observed, for example, 
damping of collective excitations [8 and inhibition of ex- 
pansion [9l [TOl [11, due to disorder. Another experiment 
[12] has examined a BEC in an incommensurate {quasi- 
random) optical lattice, in order to investigate a possible 
"Bose-glass" phase [13 . Experiments with disordered 
atomic quantum gases may provide unique insights into 
disordered quantum systems, and may uncover a rich va- 
riety of quantum phases [14 . 

Here we report experiments on a BEC of interacting 
''Li atoms subject to a well-controlled disordered poten- 
tial. While we corroborate previous transport measure- 



ments [HI [H Uni [II] , we have also probed the ground state 
density distribution and phase coherence of the disor- 
dered BEC by performing both in-situ and time-of-flight 
(TOE) imaging. While disorder inhibits transport of 
the BEC, reproducible TOE interference patterns are ob- 
served for intermediate disorder strengths Vd, reflecting 
an underlying phase coherence in the disordered BEC. 
At stronger Vd, the interference contrast diminishes as 
the BEC fragments into a "granular" condensate, which 
is expected to have no phase coherence. 

We create an optically trapped BEC of ^Li atoms in 
the mj = —1/2, m/ =3/2 state at a bias magnetic field B 
c:^ 720 G, using procedures essentially similar to those de- 
scribed in [15 . The bias field tunes the scattering length 
as via a Eeshbach resonance, to a large positive value 
(~200 ao, where ao is the Bohr radius), enabling efficient 
evaporative cooling [IH [H]. The optical dipole trap is 
produced by focusing a 1030 nm Yb:YAG laser beam to 
a Gaussian 1/e^ intensity radius of 34 /im. In the follow- 
ing, we label the axial direction (along which the optical 
trap beam propagates) as z, and the two radial direc- 
tions as X and y. At the final trap depth, the radial trap 
frequencies cOx = — 27rxl80 Hz, and the axial trap fre- 
quency iJz — 27rx3.6 Hz. Part of the axial confinement 
is contributed by a residual curvature from the bias mag- 
netic field. At this trap depth, the BEC contains ~5xl0^ 
atoms with no discernible thermal component. The cor- 
responding chemical potential /i ~ /ixl kHz, where h is 
Planck's constant. 

To create the disordered potential, we pass another 
beam, derived from the same laser as used for the opti- 
cal trap, through a diffusive plate. This results in laser 
speckle [IT] similar to that used in [SI [lO]. This beam 
propagates along the x direction. Because of the optical 
dipole force, the speckle's spatially varying intensity pat- 
tern gives rise to a disordered potential V{x,y^z) for the 
atoms. Figure [TJa) shows a 2D cross-section, V(0,y,z), 
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FIG. 1: (color) Typical ID disordered potential V{z). (a) 
False-color surface plot of the disordered speckle potential 
V{y,z). The pixel size corresponds to 2.2 /xm in the plane 
of the atoms, (b) Autocorrelation function of V{z). The line 
is a fit to e~^^^ , within a scaling factor. As expected 
for fully developed speckle, the intensity follows a negative 
exponential probability distribution. 



as imaged by a camera, while the main panel shows an 
axial cut, V{z) = 1/(0, 0,z), of the disordered potential 
through the center [x = y = 0) oi the atomic cloud. The 
disorder strength, Vd, is proportional to the intensity 
of the speckle laser, and is defined as the standard de- 
viation of the mean of V{z). The disorder correlation 
length along the z direction, a^, is defined by fitting the 
autocorrelation function {V {z)V {z Az)) to e~^^^ /"^^ 
(Fig.jljb)); we find 15 /im, which is much greater 

than the condensate healing length of ~1 /am. Since ax 
^ cFy = (Jz [IQi ilTj and the radial Thomas- Fermi radius 
of the BEG is less than 10 /im, to a good approximation 
the BEG is subject to a ID disordered potential V{z) = 
V(0, 0, z) along its axial direction. 

After producing a pure (to within our resolution) BEG, 
the disordered potential is ramped on in ~100 ms [18] to a 
specific value of Vd, which we express in terms of fi. The 
BEG is then held in the disordered and trap potentials 
while various experiments are performed, as described 
below. The final state is probed by absorption imaging 
(also along the x direction) of the atomic cloud released 
from the disordered and optical trap potentials, or by 
in-situ phase contrast imaging [T9^. 

We first describe data characterizing the transport 
properties of the disordered BEG. In one type of mea- 
surement, we slowly ramp on a magnetic field gradient 
along the z direction, which offsets the center of the har- 
monic trap from z = io z = —d. The disordered poten- 
tial is kept stationary during the trap offset. As shown 
in Figs. [2|a) and (b), without disorder the center of the 
atomic cloud follows the trap. For intermediate disorder 
strength the cloud lags behind the new trap center and is 
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FIG. 2: (color) Effect of disorder on slow transport of the 
BEG. The trap center is offset to the left for a distance d over 
a time interval of 700 ms. (a) Initial state of the BEG, without 
disorder or offset. The center of this initial cloud defines the 
reference position "0" . (b-d) Result of offsetting the trap by d 
^ 1400 /im, with the cloud subjected to a disordered potential 
of various Vd'- (b) Vd — 0; (c) Vd ~ 0.3 /x; (d) Vd ~ Images 
(a-d) are obtained by absorption imaging taken with a 1 ms 
TOF after the cloud is abruptly released from the disordered 
and optical trap potentials. Main panel shows the center by 
weight of the cloud as a function of Vd , for several offsets d. 
Vd and ji have uncertainties of '^20%. 



stretched, as shown in Fig.[2|c). For stronger disorder, as 
shown in Fig. |2|d) , the cloud is pinned at its initial posi- 
tion and does not respond to the offset. We plot the cloud 
center as a function of Vd for 4 different offset distances 
d in Fig. [2] The cloud is pinned and does not respond to 
trap offset for Vd ^ 0.8/i, within our experimental time 
scale and uncertainty. 

We have also studied the effect of disorder on the dipole 
oscillation of the trapped BEG (Fig. |3|. In contrast to 
the slow offset experiment discussed above, here the har- 
monic trap center is abruptly (2 ms) offset, and the cloud 
evolves in the presence of the shifted trap and station- 
ary disorder for a variable time before being released for 
imaging. Without disorder, the condensate undergoes 
undamped dipole oscillations in the trap with frequency 
uoz' Even disorder as small as Vd ~ 0.1 /i, however, 
damps these oscillations. Increasing the disorder strength 
to Vd ^ 0.4 /i causes the dipole motion to be overdamped, 
until when Vd ^ A^, the cloud becomes pinned at its ini- 
tial position to within our experimental resolution. 

In addition to transport measurements, we have per- 
formed in-situ [20 and TOF imaging on the disordered 
BEG, yielding information on both its density and phase. 
Figures [Ij^a) and (b) show in-situ and corresponding 
TOF images taken 8 ms after releasing the cloud. For 
intermediate Vd (Figs. |4jb) and (h)) striking random 
fringes, which we interpret as matter wave interference. 
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FIG. 3: (color) Effect of disorder on dipole excitation of the 
condensate. Data shows the cloud center as a function of time 
following a sudden axial offset of the trap, for various disorder 
strengths Vd (shown in two separate panels for clarity). The 
lines are fits to damped dipole oscillations. 



develop after sufficiently long TOF expansion of the 
BEG following release from the optical potentials. We 
note that at such Vo, the corresponding in-situ images 
(Figs. [4]^ a) and (g)), are consistent with the cloud still be- 
ing connected. By increasing Vd above /i, the disordered 
BEG becomes fragmented shown in the in-situ profile, 
Fig. [4]^i), and the fringe contrast observed in TOF di- 
minishes, as shown in Fig. [4]^j). To quantify the depen- 
dence of the TOF fringe pattern on Vd, we plot both 
the TOF and in-situ contrast vs. Vd/m Fig. [s] The 
TOF fringe contrast peaks at Vd — 0.5 /i, while the con- 
trast in the in-situ images increases up to the highest Vd 
investigated. 

Similar, but less well resolved fringes in TOF im- 
ages have been reported previously [H [11]. Our highly 
elongated BEG facilitates a systematic and quantitative 
study of the fringes. We find the positions of these irreg- 
ularly spaced fringes are reproducible in repeated mea- 
surements, even with different holding times of the BEG 
in the disordered potential. During preparation of this 
manuscript, we became aware of a closely related paper 
[2T] . which also reports reproducible fringes in TOF im- 
ages of a disordered BEG. The reproducibility suggests 
that these fringes are unlikely due to some initial (before 
release) phase fluctuations [22 in the disordered BEG, 
or due to interference of a few BEGs separated by high 
potential barriers with no well defined relative phase be- 
tween them, as in either of these scenarios the fringe po- 
sitions are not expected to be reproducible. 

Numerical simulations based on the Gross-Pitaevskii 
equation |TT1 [21] have shown that such interference can 
occur when moderate disorder-induced density fluctua- 
tions cause different parts of a phase- coherent BEG to 
expand with different velocities and overlap. As Vd is 
increased from 0, the contrast of the interference is ex- 
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FIG. 4: (color) In-situ profiles and matter wave interference, 
(a) In-situ image for Vd ~ 0.5 /i; (b) TOF image for Vd ~ 
0.5 /i. The TOF 2-D image uses a different color scale from 
the one used for the in-situ 2-D image, and the displayed 
aspect ratio of the in-situ images is reduced by 50% in or- 
der to enhance the transverse size for visibility, (c-j) In-situ 
column-density profiles on left and corresponding TOF images 
on right; (c, d) Vd = 0; (e, f) Vd ^ 0.3 /x; (g, h) Vd ^ 0.5 
/i; (i, j) Vd ^ 1.0 /i. Solid red lines are fits to Thomas-Fermi 
distributions, cuz — 27rx2.8 Hz for data shown in Figs.^and 
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FIG. 5: (color online) Contrast of the matter wave interfer- 
ence pattern in 8 ms TOF images (filled diamonds) and con- 
trast of in-situ density fluctuations (open circles) as a function 
of Vd/m- The contrast is extracted from the central axial cut 
of the column density, as its mean deviation from the Thomas- 
Fermi fit normalized by the peak value of the fit. 



pected to increase, consistent with the data shown in 
Fig. [5] at small Vd- When Vd becomes sufficiently large, 
however, the BEG fragments into multiple pieces, as is 
seen in the in-situ image, Fig.|4ji). An array of randomly 
spaced condensates should not produce visible interfer- 
ence in TOF [23l [24] , again consistent with our observa- 
tion of diminishing fringe contrast at high Vd (Figs.|4jj) 
and |5|. 
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At high Vd, we observe completely inhibited transport, 
indicating a transition to an insulator. This was seen 
in previous transport experiments, where the nature of 
the insulator is inferred to be that of a fragmented BEC 
[HI m [To]. Our in-situ images give the first direct ob- 
servation of fragmentation of a highly disordered BEC. 
Because of the exponentially suppressed Josephson tun- 
neling between fragments due to high potential barriers, 
the fragmented BEC is expected to have large phase fluc- 
tuations and no phase coherence [25J. This regime is 
analogous to a granular superconductor [3l [4] , which is 
an insulator without global phase coherence. 

Dipole oscillations are damped at moderate Vd, as 
shown in Fig. [3j and in agreement with previous observa- 
tions [8 . Excitations such as solitons, vortices and phase 
slips have been suggested [27] as possible mechanisms 
giving rise to damping. Our in-situ as well as TOF data 
at moderate Vd suggest a disordered BEC that is still 
connected (not yet fragmented) and phase coherent, in 
order to give a reproducible interference pattern. Such 
a state may provide insight into the physics of "homoge- 
nous" disordered superconducting fllms, where there has 
been tremendous interest regarding the form of the order 
parameter, possibilities of metallic states, and the fate of 
phase coherence and phase stiffness [28] during the SIT 
[3l[2^. Damped dipole oscillations were also recently stud- 
ied for a BEC in a periodic potential [29 in the context 
of a possible "Bose metal" phase [Tj. 

Strong repulsive interactions are believed to preclude 
Anderson localization (AL) [lll[25] in a BEC such as ours 
and those in previous experiments where the condensate 
healing length is small in comparison to [8l|9l HO] [11] . 
Interactions, however, can be signiflcantly reduced in ^Li 
due to an extremely gradual zero-crossing in the B depen- 
dence of as [ISl [16] , creating more favorable conditions 
to observe AL. Such experiments are currently underway. 
A similar zero-crossing has been identifled in ^^K [26 . 
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